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Abstract 

Although the statistical mechanics of periodically driven ("Floquet") systems in contact with a 
heat bath has some formal analogy with the traditional statistical mechanics of undriven systems, 
closer examination reveals radical differences. In Floquet systems all quasienergies ej can be placed 
in a finite frequency interval < Sj < uj (with uj the driving frequency and h = 1), and the number 
of near degeneracies {{ej — £k\ < S for arbitrarily small 6) in this interval grows without limit as 
the dimension of the Hilbert space increases. As we noted in a previous paper, this leads to 
pathologies, including drastic changes in the Floquet states, as A^ increases. In earlier work on 
Floquet systems in contact with a heat bath these difficulties were put aside by fixing A^, while 
taking the coupling to the bath to be smaller than any quasienergy difference. This led to a simple 
explicit theory for the reduced density matrix, but with some major differences from the usual time 
independent statistical mechanics. We show that, for weak but finite coupling between system and 
heat bath, the accuracy of a calculation within the truncated Hilbert space spanned by the N 
lowest energy eigenstates of the undriven system is limited, as A^ increases indefinitely, only by the 
usual neglect of bath memory effects within the Born and Markov approximations. As we seek 
higher accuracy by increasing A^, we inevitably encounter quasienergy differences smaller than 
the system-bath coupling. We therefore derive the steady state reduced density matrix without 
restriction on the size of quasienergy splittings. In general, it is no longer diagonal in the Floquet 
states. We analyze, in particular, the behavior near a weakly avoided crossing, where quasienergy 
near degeneracies routinely appear. The explicit form of our results for the denisty matrix gives a 
consistent prescription for the statistical mechanics for many periodically driven systems with N 
infinite, in spite of the Floquet state pathologies. 
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I. INTRODUCTION 



Time-periodic quantum systems are of great importance in several fields of science: 
atomic and solid-state systems driven by monochromatic electromagnetic fields quan- 
tum ratchets 2]; quantum chaos j^. The Floquet theorem gives the solutions of the time- 
dependent Schrodinger equation with a time-periodic Hamiltonian (period l-njuj) in the 
form exp(— with uii) exhibiting the same periodicity as the Hamiltonian. The 
"quasienergy" e plays a role analogous to the energy in static systems, determining in this 
case the behavior of the wave function under finite one-period, rather than infinitesimal, 
time translations. The quasienergy e is defined only mod cj, and we follow the convention 
of taking < e < u;. 

Real Floquet systems will ordinarily be in thermal contact with a bath. Energy is being 
continuously supplied by the periodic driving field, with relaxation provided by the bath. 
This leads to a fundamental question of statistical physics: Is there an asymptotic, time 
periodic, ensemble characterizing the quantum system and, if so, how does it depend on the 
coupling to and the temperature of the bath? For a time independent system with Hamilto- 
nian Hq the answer is given by the canonical reduced density operator p = exp{—f3HQ)/Z, 
where f3 is the inverse temperature, and Z is the partition function. This density operator is 
independent of the form and magnitude of the coupling to the environment, to lowest non- 
trivial order in that coupling, and the operator is diagonal in the energy eigenstates of the 
system Hamiltonian. Thus, even if two states are separated by an energy small compared 
to their coupling with the bath ("near degeneracy"), or even if they are degenerate, their 
statistical weights are given by their energy, and there is no anomaly associated with their 
near or strict degeneracy. That simplicity does not carry over to the case of near degeneracy 
in quasienergy . 

For time-periodic systems in thermal contact with a bath the general long time behavior 
has not yet been fully clarified. In particular, the inevitable appearance of near degeneracies 
of quasienergies has not been addressed. Earlier work includes a number of papers 0, 

which have discussed specific physical systems, treating either those whose quantum 
states span only a finite space, or others whose Hilbert space has been truncated to a finite 
set and then studied numerically, but always with coupling to the bath small with respect to 
any quasienergy spacing. Within such a finite state space, and in the limit of weak coupling 
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to the bath, a description of driven dissipative systems, as apphed to driven Rydberg atoms 
in thermal contact with discrete wave guide modes, was introduced by Bliimel et al. j4|, 
and extended by Kohler et al. % to driven one-dimensional harmonic and anharmonic 
oscillators. An extension to strong coupling to the bath was presented by Thorwart et 
al. 1^ for tunneling and relaxation of a particle in a bistable one-dimensional well, but 
using methods and approximations suitable only when a (small) finite number of system 
energy eigenstates is significantly accessible. A discussion of the statistical mechanics of a 
more general class of periodically driven systems was given by one of us ll|], restricted to 
a finite Hilbert space with a non-degenerate quasienergy spectrum, in the limit of coupling 
to the bath much weaker than the smallest quasienergy spacing. In this limit, where the 
density operator is assumed diagonal in the representation of the Floquet states, the matrix 
elements can be obtained by solving a set of linear rate equations. In sharp contrast to 
the time independent case, they do depend on the explicit form (but not the magnitude) 
of the coupling between system and bath, and they are not simple smooth functions of the 
quasienergy. 

In the general case of time-periodic quantum systems with a Hilbert space whose num- 
ber of dimensions increases without limit, the situation is more complex. Since the 
quasienergies are restricted to a finite range of width tu, the mean spacing between neigh- 
boring quasienergy levels approaches zero as their number increases. Therefore, any fixed 
coupling to the thermal environment is never small relative to all level spacings for a suf- 
ficiently large basis size. It becomes essential to deal with degenerate and near degenerate 
quasienergies, a situation which has been explicitly excluded in previous analyses. Indeed, 
in the limit A^ — cx) an infinite number of quasienergy degeneracies avoided by arbitrar- 
ily weak level repulsions may occur, and many anomalies of the spectrum and the Floquet 
states appear, as was pointed out in [l2|. So it is of special importance to determine the 
effect of the bath when the coupling becomes comparable to or even large with respect to 
one or more of the relevant quasienergy splittings. 

In this paper we address the critical question of the appropriate description of the steady 
state of time-periodic quantum systems with small but finite coupling to the bath. We 
review in Sec. II the derivation from the master equation of the set of linear equations that 
determine the reduced density matrix. This matrix is in general no longer diagonal in the 
basis of the Floquet states of the system, but it does approach a limiting time periodic form. 
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We obtain a set of linear homogeneous equations which determine all of the density matrix 
elements in the long time, stationary time-periodic state, Eq. (12.201) . This result is no longer 
limited to the regime of bath coupling weak compared to all quasienergy spacings. For a 
finite Hilbert space and in the limit of infinitesimal coupling we discuss the differences from 
the conventional statistical mechanics of time-independent systems. 

We next study in detail the influence of a single weakly avoided crossing of Floquet 
quasienergies (as a function of driving field strength) on a finite density matrix, since it is 
there that the near degeneracies of interest develop. In Sec. Ill we find that the important 
consequences of such an avoided crossing can be incorporated by introducing an additional 
effective coupling (13. 4p between the two states involved, mediated by the bath. This coupling 
decreases with distance from the avoided crossing (which conveniently parameterizes the 
quasienergy separation of the two states), and it depends on the ratio of the coupling strength 
r of the two states to the bath, and their minimum quasienergy splitting, A. Remarkably, 
all of the diagonal elements of the density matrix in the representation of the "diabatic" 
states defined well away from the center of the avoided crossing, a basis which is fixed 
throughout the crossing, are determined by a set of linear rate equations (13.31) of exactly 
the structure familiar from ordinary time-independent statistical mechanics, but with this 
single added effective coupling rate. For F ^ A the additional coupling influences not 
only the (diagonal and off-diagonal) reduced density matrix elements for the two states 
themselves, but may completely change the whole steady state reduced density matrix, so 
that average system observables will in general change radically with driving field strength 
in the range of the avoided crossing. Physically, high order near resonances of the driving 
field between eigenstates of the undriven system can populate high energy levels, which 
then affect populations of lower energy levels in the relaxation cascade. In the opposite 
limit, F 3> A, the coherence necessary for the high order resonance is disrupted by the 
interactions with the bath, and the avoided crossing has negligible influence on the reduced 
density matrix which, in this limit, is approximately diagonal in the diabatic basis. We 
present a simple, albeit somewhat artificial, example which can be solved in some detail, 
specifically exhibiting these effects. 

We then use this insight to discuss time-periodic quantum systems with a Hilbert space 
of increasingly large dimension in Sec. IV. We find that weakly enough avoided crossings, 
while having great impact on the Floquet states, have no significant consequences for the 
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density matrix (in a A-independent basis). It is this central result which justifies the neglect 
in many cases of all but a finite basis set in the calculation of the statistical mechanics of 
a Floquet system, as is essential for any numerical study. Specifically, for the calculation 
of statistical properties, we establish the validity of using the conventional approximation 
of eliminating from the undriven system spectrum all states of sufficiently high energy (the 
cutoff to be determined by the accuracy desired), and allowing the periodic driving field to 
mix only the remaining states. Within that finite basis set we have derived a prescription 
for calculating the full reduced density matrix, in general no longer diagonal in the Floquet 
representation. A brief summary is given in the conclusion. Sec. V. 



II. REDUCED DENSITY MATRIX 



A. Equations of motion and time periodic solutions 

We review and extend the derivation of the steady state solution for the reduced density 
n,at... in the waak-coupHng Unnt ( "Bo.n-Mark„v app.oxi„.atio,0 fo. F,o,uet .y=te,n=, 

previously discussed 17|] only for situations where all quasienergy spacings were large 
compared with the coupling to the bath (and therefore only for systems restricted to a finite 
number of states). In contrast, we stress that the validity of the more general equations is 
not limited to that regime, and we give a detailed justification for an approximation using 
time-averaged rates. 

We write the Hamiltonian generally as the sum of Hs{t), describing the dynamics of the 
system driven by a time-periodic external field with period 2tt/uj, the Hamiltonian Hf, of the 
bath, and the interaction between system and bath Hgh, 

H{t) = H,{t) + H, + H,,. (2.1) 

The total density operator of this closed coupled system, prit), is the solution of the equation 
of motion, 

tpT{t) = [H{t),pT{t)]. (2.2) 
We trace over bath variables to obtain the reduced density operator for the system, 

p = Tib pT, (2.3) 
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and do the standard recasting of the master equation into an integro-differential equation 
(see, e.g., Q, Q), whose iteration gives directly the expansion in powers of the system-bath 



interaction: 

d /■* 
d^' 



~p{t) = - f dt'Ti^ \H,,{t), \Hs,{t')rPT{t')\\ . (2.4) 

Jo L L J J 



Here the tilde labels operators in the interaction picture. We have omitted the linear term 
involving Trb-ff^fc, which ordinarily vanishes. This is the case, e.g., for an interaction between 
system and bath characterized by creation or destruction of single bath bosons, such as 
phonons or photons, or of free electrons or other fermions. Otherwise, this equation is still 
general and exact. 

In order to make further progress, which depends on being able to trace over the bath 
variables under the integral, we make two essential simplifications: (i) second order pertur- 
bation theory [igI in Hsb and (ii) the Markov approximation. Their validity relies ultimately 
on being able to neglect any memory effects within the bath. The correlation times associ- 
ated with the bath must be small compared with the times characterizing the evolution of 
the reduced density operator (see the condition (12.101) below). First, to lowest non- vanishing 
order in ifsf, we take the total density operator prii') under the integral as a simple product 
of system and bath factors, 

Prit') ^ p{t')p,. (2.5) 

The physical meaning of this approximation is that we neglect more than a single interaction 
of the system with the bath during the bath correlation time. The effect of higher order 
correlations, where subsequent interactions depend on the state of the bath modified by 
earlier ones, is negligible, as we show in more detail in Appendix \M 
We make the specific, but inessential, choice of couphng Hamiltonian: 

Hsi, = lAB, (2.6) 

where the Hermitian operator A acts only on the Floquet system, and the Hermitian operator 
B acts only on the heat bath. A sum of such factorized operator products, 'jJ^aAaBa, is 
also readily treated. But if the various bath operators are dynamically uncorrelated 
with each other, then the results below are modified only by an extra index on each matrix 
element of A and each spectral function g{E) (to be defined below), and a sum over that 
index. We take the single product to avoid the consequent algebraic clutter, since no essential 



7 



difference results in the ultimate equations (the extra sums being merely incorporated into 
the definition of the complex "rates" below (12.141) ). Then the trace over bath variables 
involves only the single combination 

G(r) = Trfcp,5(t)S(t-r), (2.7) 

This bath correlation function is characterized by a decay time Tc, which we assume to be 
much shorter than the times characterizing the evolution of the reduced density operator 
p{t) (see Eq. (12.101) below). Then the main contribution to the integral in (12. 4p is from values 
of t' within a range of order Tc below the upper limit t. Thus to a good approximation we can 
extend the lower limit to — oo and also make the Markov approximation: replace t' by t in 
the argument of pr under the integral. Returning to the Schrodinger picture and changing 
the integration variable to r = t — t', we then have 

^ + imt),pit)] = -7' dr {G(r) [AAit - r, t)p(t) - i(t - r, t)pit)A] + H.c.} . (2.8) 

The rate of change of the reduced density operator due to the coupling to the bath is of 
order 

T,=j\{A'){B'), (2.9) 
suggesting that the condition for validity of the "Born" and Markov approximations is 

TeFs < 1. (2.10) 

Appendix |X] demonstrates that this same inequality assures that higher order corrections 
are negligible. The validity of the above approximations is discussed in the literature on 



undriven systems; a particularly clear discussion is given in Ref. [15|] . As the approximations 
have been made within the interaction picture, they are equally valid for a time-periodic 
system Hamiltonian Hs{t). We emphasize that we do not require infinitesimal coupling to 
the bath. That coupling need only be small enough to satisfy the condition (I2.10p . 

We recall that the solutions of the isolated system Schrodinger equation idt\(p) = Hs(t)\(f)) 
are of the Floquet form \(pj(t)) = exp{—iejt)\uj(t)) with the quasienergies ej chosen to lie 
in the strip < e < uj and with the time-periodic parts \uj{t)) = \uj{t + 211/00)) forming 
a complete orthonormal basis for the system at any given time t. We introduce the matrix 
elements of the reduced density operator in this time-periodic basis, 

A,(t) = (w,(t)|p(t)|«,(t)), (2.11) 



which will turn out to be very convenient in the following calculations and which, as we will 
show later, leads to time-independent density matrix elements in steady state. In this basis 
the equation of motion of the reduced density matrix elements takes the form 

(^ + ^^^.)a.W = -7' E e^(^+"^)-*[ + A,(M)A,Km)py(t) / drG(r)e^(^--'"i^.l2) 

'^^ k,l,M,m ''^ 

/•oo 

- A,k{m)pM{t)MM) / drG(r)e*(^--™-)^ 

Jo 

+ p,k{t)Aki{m)MM) / drG'*(r)e*(^--'"-)- 

^ 

where we defined the Fourier components of the matrix elements of the system operator A, 

{u,{t)\A\u,{t)) = Y^A,,{m)e'''^\ (2.13) 

m 

with the property Aki{m) = A^i.{—m), and where we have introduced the notation Eij = 

Si Sj. 

We simplify the notation by introducing complex time-periodic "rates" 

R^rMit) = J2 e*'^"*i?.j;fc/(i^), (2.14) 

K 

with Fourier coefficients 

R^rMiK) = Y.R'l];kiiK) = 2n^'Y.M^ + K)Al,im)gieik-rnuj), (2.15) 

m m 

where g{E) is the Fourier transform of G{t). We will approximate 

dTG{T)e-'^^ ^7Tg{E), (2.16) 

neglecting a principal value contribution. From the definition of G{t), Eq. (12. 7p . one finds 

g{E) = J2Pu\{'y\B\iy')\^6{E + E,,-E,), (2.17) 

where u labels bath states of energy E^ which are in thermal equilibrium and have occupation 
oc exp{-(3E^). 

Then we have as the fundamental equation of motion of the reduced density matrix 
elements in the Floquet representation, the following set of linear differential equations: 

{^^+^^v)p^At) = -lll[piAi)RikAt)+Pum*,k;ikit)-PH^ (2.18) 



This is the generahzed Master equation for all density matrix elements. We next replace 
the time periodic rates Rij;ki{t) by their average over one driving period, 

RirM = E R'l^-M = Rir,ki{K = 0). (2.19) 

m 

I — I 

This approximation [3] implicitly assumes that the density matrix elements do not vary 
substantially over a period of the driving, and it is further discussed in Appendix |Bl We 
then find a steady state solution of ( 12.18^ for which the matrix elements pij are independent 
of time: 



- E [PljR'ik;lk + PilR*jk-lk - Pkl{Rlj;ki + Rki;lj) ■ {2.20) 



i^ijPij „ 

^ k,l 

B. Large quasienergy splittings 

The states of some systems of interest span only a finite dimensional space (e.g., spin 
systems). Because of the finite number of states, at a given driving strength A there will 
be a smallest quasienergy separation. In this subsection we take that smallest separation 
to be large compared to the coupling to the heat bath (an assumption made, in fact, in all 
earlier work on driven dissipative systems in the weak-coupling limit), and we will point to 
both the similarities to and the differences from the conventional statistical mechanics of 
time-independent systems. 

In the limit where the absolute values of all "rates" \Rij-ki\ are sufficiently small compared 
with all splittings \Bij\, it is clear from f l2.20p that all off-diagonal density matrix elements 
can be neglected, and we arrive at the familiar form, 

= Pii'^Rik — PkkRki, (2-21) 
k k 

relating the state occupation probabilities Pi = pa. Here the diagonal time-independent 
rates, 

Rki = Rki;ki = 2TT-f'^'^\Aki{m)\^g{eik - mu), (2.22) 

m 

are of the form of the usual transition rates of perturbation theory, the square of coupling 
matrix elements times the density of states (at the energy transfered from bath to system), 
summed over initial states of the bath with appropriate thermal equilibrium weighting. Of 
course, in the static case the equation also holds without the summation over k (detailed 
balance). That is not the case here. 
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As was pointed out earlier [o], [ll|, there are some important differences in this result 
from that familiar in static problems. Because of the possibility of absorption and emission 
of multiple quanta of the time- dependent driving field, the occupation probabilities are not 
simple Boltzmann factors, and indeed they depend on the detailed form of the Hamiltonian 
coupling the system to the heat bath. From the definition 02.1 7p of the spectrum of the bath 
correlation function we have 

g{E) = g{-E)e-P^, (2.23) 

which is used in the time-independent case to derive the detailed balance equation. In a 
Floquet system it leads, using fl2.17p and (12.221) . to the relation 

i?™, = i?-™e^("^-"^'+'""), (2.24) 

relating only the pieces of the transition rates between system states j and j' which involve 

the transfer of exactly m quanta from the field. After summation over m we find for the 

ratio of the full transition rates between these states, 

p pm 

^n' ^ ^™ jj' ('2 25) 

There are situations where the i?^, have for some pairs (j, j') just one value of m, say mo, 
with a non-negligible value. This is the case for a very small driving amplitude (A ^ 1) or 
for high-lying states of the undriven system that are weakly affected by the periodic driving, 
if the levels j and j' are far enough from any avoided crossing. Then in the above ratio the 
terms i?^? cancel and the form of the coupling Hamiltonian if^fc is no longer relevant. One 
gets Rjj'/Rj'j ~ exp[P{ej —Sj' + mQUj)], much as in the time independent case, but reflecting 
the need for itlq quanta of the driving field. 

In general, the states are mixed by the driving field, giving rates R'j]' where more 
than one value of m contributes significantly. Consequently the above ratio Rjj'/Rj'j will 
no longer depend just on a single quasienergy difference. In order to find the stationary 
weights pj = pjj one always has to solve Eq. (12.211) . This solution will in general not fulfill 
the detailed balance equations. Even for the simple situation above, where the direct and 
inverse rates for states j and j' have a simple Boltzmann factor ratio, we can not then 
conclude that the corresponding probabilities pj and pj' have the same ratio. As in the case 
of a time independent Hamiltonian, the stationary solution of Eq. (12.211) will not depend 
on the overall amplitude 7 of the coupling Hamiltonian (if 7 is sufficiently small for lowest 
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order perturbation theory to be valid). In contrast to the time-independent case, however, 
the solution will depend on the rates RJji and thus on the precise form of the coupling 
Hamiltonian Hgh, as noted in 



III. ISOLATED WEAKLY AVOIDED CROSSING 

As a function of the strength of the driving field A the quasienergy spectrum of Flo- 
quet systems is generally pervaded by weakly avoided crossings, near which quasienergy 
separations are small and the Floquet states are rapidly varying. These features can have 
important consequences for the density matrix. We consider a finite set of Floquet states 
for which, in the neighborhood of the driving field strength A being studied, there is a single 
weakly avoided crossing, with small quasienergy splitting A at the center of the crossing 
at Aq near A. At the same time all other quasienergy separations at A are assumed to be 
large compared to the appropriate measure of the coupling strength to the bath. Then only 
the two Floquet states involved in the weakly avoided crossing are sensitive functions of A 
over the small crossing region. For very weak coupling to the bath compared with A we 
know from the results of Sec. JIB that p is diagonal in the basis of those Floquet states. In 
contrast to the time independent case, where the density matrix elements depend only on 
the state energies, here it will change rapidly as the Floquet states do, not only within the 
subspace of the states of the avoided crossing but, as we will show, in general throughout 
the rest of the density matrix, as well. In the opposite limit, when the system-bath inter- 
action effectively broadens those two states by much more than A, we will show that there 
is little influence from the crossing on the density matrix (in a A-independent basis) or on 
the statistical properties of the system determined therefrom. This will allow us in the next 
section to extend our analysis to many driven systems for which the quasienergy spectrum 
is dense, as is often the case in practice. 

We review briefly the description of the states and quasienergies in the neighborhood 
of such a weakly avoided crossing, introducing the notation which will be used later. We 
assume that within a suitably small neighborhood of driving strengths around Aq the time- 
periodic parts \ui{t)) and |^^2(^)) of the unperturbed "diabatic" states (i.e., if there were no 
avoided crossing) do not depend significantly on A, and that the unperturbed quasienergies 
are locally linear in A with slopes ai and o"2, respectively. Then, the quasienergies and the 
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time-periodic part of the Floquet states can be well described by the results of diagonalizing 
a 2 X 2 Hamiltonian matrix (the Floquet analog of nearly degenerate perturbation theory, 
using the extension of Rayleigh-Schrodinger perturbation theory developed by Sambe 
The quasienergies are given by 

ea,t = eo + ^^i±^(A - Ao) ± ^ {a, - a^nx - Xo)' + A' , (3.1) 

where the indices a, b refer to the two "adiabatic" Floquet states continuous in A, and A 
is the minimum quasienergy splitting, at the center of the crossing, where A = Aq. This is 
shown schematically in Fig. 1. The adiabatic states a and b exchange their character when 
A is varied, while by assumption the diabatic states 1 and 2 remain unchanged. 
We define the dimensionless distance d from the avoided crossing at A = Aq as 

d=(A-Ao)(ai-a2)/A, (3.2) 

the ratio of the quasienergy separation of the unperturbed states at A to the gap A. The 
central region of the avoided crossing is thus given by \d\ < 1. 




FIG. 1: Behavior of quasienergies around an avoided crossing of size A. The dashed Hnes show 
diabatic (unmixed) states 1,2 and the solid curves the adiabatic Floquet states a,b. 

A. Diabatic basis and effective rate i?^^ 

Over the narrow range of the weakly avoided crossing the diabatic states \ui{t)) and 
\u2{t)), as well as the time periodic parts of all other Floquet states, are approximately 
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independent of A, so they constitute a natural fixed basis of time periodic functions to 
use in studying this region. In Appendix [C] we derive the set of rate equations for the 
reduced density matrix elements pij in this diabatic basis. Quite strikingly, under a number 
of plausible assumptions, they allow for a separate determination of the diagonal and off- 
diagonal elements. The former are simply given by 

= Pii E ~ E Pjj^ii^ (3-3) 
i j 

where the rates designated with a bar, Rij, are defined in the diabatic basis in analogy to 
Eq. (223), except for 

R,, = R,, = R- ^ 4^,^(r/A)2' ^2-^) 
with the superscript on R'^'^ a reminder that this is an effective rate associated with the small 
quasienergy separation near the avoided crossing (ac). Here we have introduced the symbol 

r = ^ {Rik + R2k) + -Rii + R22 — 2-Rii;22, (3.5) 

k 

where the prime on the summation denotes the restriction k 7^ 1,2. (Note that the sum 
of the last 3 terms is positive, in spite of the minus sign; by definition they can be written 
27r7^ X^m [^11 ('^) ~ A22{jn)\'^ g{—muj) .) The rate V will turn out to be the relevant coupling 
to the heat bath, specific to the crossing of states 1 and 2, that will be the measure of 
whether that coupling can be considered weak or strong relative to the quasienergy splitting 
of those states. 

These equations (13. 3p for the diagonal density matrix elements pij in the diabatic basis 
have exactly the formal structure of the corresponding equations in the limit of bath coupling 
small compared to all quasienergy splittings, Eq. (I2.2ip . The results differ from that simpler 
situation in three essential ways: 

(i) The effective rates R'^^ between the two special states 1 and 2 are not simply the 
second order rates imposed directly by the bath. They reflect the mixing in a more complex 
way, effectively including the driving induced coupling to all orders, and they can be simply 
determined from the rate F defined in Eq. (13.51) and the properties of the avoided crossing, 
namely the distance d and the minimal quasienergy splitting A. 

(ii) The dependence on the parameter A, which in the original adiabatic basis strongly 
affects a great number of rates, appears in the diabatic basis only in the single rate R^"^ (see 
Eq. (13.41) . where the distance from the avoided crossing, ci, is linearly related to A). 
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(iii) These are not the only density matrix elements; there exist also two non-vanishing 
off-diagonal elements pu and p2i = Pu in this diabatic basis, given by 

Pll — P22 /„ «\ 

Their importance can be assessed by a comparison of the absolute value of pu to the sum 
of the positive and real values pu and p22'- 

IP12I IP11-P22I 1 ^3^^^ 



Pll + P22 Pll + P22 J4rf2 + (r/A)2 r 



which is negligible in the limit F ^ A, (that is, where the coupling to the heat bath is 
stronger than the minimum quasienergy separation of the Floquet states). This demonstrates 
that in this limit the density matrix in this two-dimensional subspace is approximately 
diagonal in the diabatic basis. In the opposite limit, the density matrix in the subspace is 
nearly diagonal in the adiabatic Floquet basis. The appropriate basis to diagonalize p is 
given in general in Appendix [Dl 

Although we have no general analytic solution of the equations (13. 3p for the diagonal 
density matrix elements, we can say a good deal about the consequences of a weakly avoided 
crossing, as the relevance of the effective rate R'^'^ depends crucially on the ratio of the bath- 
induced rate F to the size of the avoided crossing A. 



1. r < A 

In this limit (which is just the case discussed in Sec. IIB above) the rate F characterizing 
the coupling to the heat bath is even smaller than the quasienergy splitting at the weakly 
avoided crossing. Consequently, the effective rate R^, Eq. (13. Sp . is much larger than F in 
the center of that crossing, 

R^^{d = 0) > F. (3.8) 

The immediate consequence of this large effective rate is that the diagonal density matrix 
elements pu ~ P22 become equal there, which may drastically change all other weights pjj, 
coupled directly or indirectly via Eq. (13. 3p to pu or P22, relative to the corresponding values 
just outside the avoided crossing. An explicit example demonstrating such a drastic change 
for the entire density matrix is studied in Sec. IIII B[ 
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The off-diagonal matrix element pu is not necessarily small, though, since it is given 
by the product of one large and one small factor; see Eq. (13.61) . In fact, as discussed in 
Appendix [D], the off-diagonal element is negligible in the adiabatic basis only. 

2. r » A 

In this limit the rate F induced by the coupling to the heat bath is stronger than the 
size of the weakly avoided crossing. Consequently, the effective rate R^^, Eq. (13.51) . is much 
smaller than F throughout the whole range of the avoided crossing, 

R^' < F. (3.9) 

This suggests that the effective rate R^'^ due to the avoided crossing has negligible influence 
compared to the other rates affecting states 1 and 2. Therefore the entire density matrix 
is not affected by the existence of the weakly avoided crossing. The off-diagonal matrix 
element pu is negligible in this limit. 

There is an exception to the above conclusion if the two sums X/fc -^lA; and X^fc R2k con- 
tributing to F have vastly different magnitudes. This is possible, as we are explicitly looking 
at a very weakly avoided crossing, and this typically arises from such dissimilar pairs of 
states, such as those arising from undriven states one of which is of low and the other of 
very high energy. In what follows we will take state 2 to be the one with higher average 
energy {Hg). There are two cases: 

• J2'k^2k R^"^ ^ J2'k^ik- This case may occur if the heat bath is dominated by 
phonons (or by other bosons) with an upper energy limit Ud, small compared to the 
separation in average energy of the diabatic state 2 and its neighbors, so that state 2 
will decay primarily to state 1, with flux p22R^- Similarly, state 2 is fed primarily by 
the direct process from state 1 with flux puR^'^ (in the extreme limit of all Rj2 = 0, 
only by this process). Since stationarity demands that the incoming and the outgoing 
flow for state 2 are equal, one finds pu ~ P22, with corresponding large changes in 
other state occupations pjj. This modification of the density matrix occurs for all 
values of d and F/A for which the presumed rate inequalities hold. Values of those 
parameters can be significantly larger than unity, but eventually large enough values 
will reduce R^'^ to the point where the lower inequality is violated. 
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• J2'k^ik ^ H'kR^k'- This case may occur when the heat bath is dominated by 

the electromagnetic radiation field, with photons of arbitrarily high energy and large 
density of states. Then state 2, essentially a high energy undriven state, will decay 
rapidly. By precisely the same argument used in the previous case, pn ~ p22; and 
states now strongly connected to 2 acquire substantial occupancy. Again, the ranges 
of d and F/A over which this situation holds can be large. 



B. A simple example 

It is instructive to obtain specific results for a simple, albeit somewhat artificial, example 
of a system. We assume first that for Floquet states labeled by some index, n = 1,2, ... N, 
the thermal bath leads to coupling only of neighboring states — i.e., the only non-zero rates 
are of the form Rn,n±i, as shown in Fig. 2, and that we are in the limit of coupling to the bath 
weak with respect to all quasienergy differences. In this special situation the lowest state 1 
has to be in detailed balance with state 2 and from this it follows that all states are in detailed 
balance. If we make the specific choice Rn,n-i = Rd and -R„,„+i = Ru, each independent of n, 
and denote the ratio Ru/ Rd = 7 < 1, then the occupations p„ = p„„ = (1 — 7)7"^V (1 ~7^) 
fall off exponentially with the index n. We assume to be sufficiently large that we can 
neglect (for algebraic simplicity) the factor 7^ below. 

Now we want to study how these occupations change if there is a single weakly avoided 
crossing of states rii and n2 (in the previous section these were the states called 1 and 2). 
This introduces the additional rate R^. It simplifies the algebra to define the fiuxes between 
states n and n — 1, 

Fn,n-l = PnRn,n-l — Pn-lRn-l,n, (3.10) 

in terms of which we can rewrite the basic equations derived above for the diagonal density 
matrix elements, (13.31) . as 

Fn2,n2 — 1 Fn2+l,n2 (3-11) 
-Fni+l,ni — -Fni,ni-1 (3-12) 

Fn+i,n - Fn,n-i for u ^ ni,n2 and n^l (3.13) 
F2,i. (3.14) 



(Pni -Pn2)R''^ = 

= 
= 
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Their solution is 



; n < ni 



""n.n-i = < F ; Til < n < n2 



(3.15) 



with 



; n > n2 



(3.16) 



being the flux from rii to n2 due to the avoided crossing. This result is illustrated in Fig. 2. 
The states above n2 and the states below ni are still in detailed balance — i.e. -F„,n_i = 0, 
as without the avoided crossing. The states from rti to n2, however, are no longer in detailed 
balance and there is a stationary flux F connecting them, due to the avoided crossing. 



n,+l 




n,+l 



FIG. 2: Simple example with a thermal bath coupling neighboring states only (double headed 
arrows) and flux F introduced by the effective rate due to an avoided crossing of states ni and 
n2- 

For the occupations Pn we find 



Pn 



Pni 



r+7" 



r+1 



Pn2 7 



Hi < n < n2 



n > 77-2 



(3.17) 
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with 



and pi determined by the normahzation J2nPn = 1- 




FIG. 3: Dramatic change in occupations pn with (sohd curve) and without (dashed curve) avoided 
crossing. 

The parameter r measures the ratio of the effective rate i?^"^ describing the direct con- 
nection of the state n2 to the state with which it is crossing, namely rii, relative to its rate 
of loss Rd to the state directly below it, times a factor of order one. There are three distinct 
interesting ranges of r in terms of the behavior of the density matrix elements: 

(i) When r ^ 1 we have p„2 ~ Pm and almost constant weights between rii and n2. The 
weights of states above n2 are much larger than they would be without the avoided crossing, 
due to the drastically increased weight Pn2, which is being fed directly by rii. 

(ii) For ^"2^"i ^ r ^ 1 we have Pn^ ~ Pmf, still a much higher occupancy than in the 
absence of the avoided crossing. Again the weights of states above n2 are much larger than 
they would be without the avoided crossing. This case, with r only slightly less than 1, is 
sketched in Fig. 3. 

(iii) If r decreases still further, r <^ ^'^-a-ni^ then p„ ^ pi'j^~^, which is the result without 
the avoided crossing. Equilibration with the bath will be completely dominated by the 
transition processes along the ladder of states. 

This model admittedly has some special features which make it soluble. For instance, 
there is detailed balance in the absence of an avoided crossing. But it illustrates the dramatic 
change in all occupations due to a single avoided crossing. 
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IV. LIMIT OF INFINITE DIMENSIONAL HILBERT SPACE 



As we have already noted, there are some interesting Floquet systems, such as finite spin 
systems, which are hmited to a space spanned by a finite number of basis states. Their 
statistical mechanics can then be calculated (numerically) from the equations and analysis 
of the previous sections of this paper (and we will not be concerned with them further here). 
However, many others, like a particle confined by a potential well, have eigenstates spanning 
infinite Hilbert spaces, with a quasienergy spectrum which becomes dense and Floquet states 
which change rapidly as the basis size N increases, so that the above results do not readily 
apply. In fact, our results also permit the calculation to any desired accuracy of many of 
these, as well, by appropriate truncation of the basis set, as we now show. 

We limit ourselves to systems with an unbounded discrete spectrum of undriven energy 
eigenstates. We take the truncated Hilbert space spanned by the N lowest energy eigenstates, 
the usual choice made in numerical calculations. Within this approximation the Floquet 
states of the isolated system (at any fixed time) are confined to that finite dimensional space. 
For the undriven system the statistical (Boltzmann) weights of the states neglected by the 
truncation can be made arbitrarily small by a suitably large choice of the energy cutoff. We 
now increase the Hilbert space accessible to the driven system by the addition of the (A^+l)st 
undriven energy eigenstate. With increasing A, as its quasienergy nears that of one of the 
original Floquet states, it will mix significantly with that state. But we have shown above 
(section [HI A 21) that, if the corresponding avoided crossing is sufficiently weak, A ^ F, then 
the reduced density matrix is approximately diagonal in the unmixed diabatic states, and the 
addition of the new state to the basis set remains of little statistical significance. Therefore, 
if we can demonstrate that all of the states neglected by the truncation of high energy 
undriven eigenstates mix with those retained only through such weakly avoided crossings, 
then we can choose the energy cutoff of the basis to obtain statistical properties to any 
desired accuracy. 

Indeed, this desired behavior of the quasienergy of a Floquet state arising from a high 
energy undriven eigenstate has been seen in numerical calculations Such a state is 
found to mix with the Floquet states that arise from low energy undriven eigenstates, as 
A increases, only near weakly avoided crossings. In the appendix of Ref. 12| we have 



analytically demonstrated the effect explicitly, with an exponential decrease with basis size 
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of newly introduced quasienergy gaps, for systems whose undriven energy eigenvalues 
above some level exhibit increasing spacings between successive levels (such as a particle 
confined by a one- dimensional potential which increases faster than harmonically). This 
can be rationalized in terms of a perturbation expansion: an increasingly large number of 
driving field quanta is needed to connect to higher lying states, which are modified only in 
increasingly higher order perturbation theory. More generally, numerical studies suggest 
that once the first undriven states have been included in the basis, the Floquet states 
arising from the first A^' of these, for suitably large N — N', are affected by further extension 
of the basis set only by weakly avoided crossings. Again, perturbation theory suggests that 
this is a reflection of the very different spatial behavior of these states. 
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For a class of periodically driven systems in an infinite Hilbert space we earlier found 
many anomalies due to the infinite number of ever more weakly avoided crossings. The 
present results show that these anomalies, and the non-convergence of the Floquet states 
with increasing basis size, lead to negligible consequences for the reduced density operator, 
or for the statistical properties calculated therefrom, when the system is in fixed weak contact 
with a thermal bath. 

In practice, one has to determine the quasienergy spectrum for increasing basis size A^ 
until (i) newly appearing avoided crossings are small enough, A ^ F, and (ii) the statistical 
(Boltzmann) weight of the newly introduced states, plus that of all neglected states, is small 
enough to achieve the desired accuracy. Once the basis size A^ is appropriately chosen, 
one can use Eq. fl2.20p to determine the density matrix elements and calculate statistical 
properties. 

It is important to recognize that, because of the requirement A ^ F for avoided crossings 
with all neglected basis states, the minimum truncation size A^ depends in general on the 
magnitude, as well as the form, of the system-bath coupling, since both enter the deter- 
mination of F. Indeed, the full statistical properties of the Floquet system will depend on 
that magnitude (we have already emphasized in Sec. IIB the general dependence on the 
form of the coupling), even if the coupling is very weak: For a given magnitude we choose 
A^ to achieve a desired accuracy. We now reduce the magnitude (keeping the form fixed) 
sufficiently that F becomes less than one of the gaps A associated with an avoided cross- 
ing of a state neglected by the truncation with a state that has been retained. Then that 
avoided crossing becomes no longer irrelevant. It leads to a large rate R^'^, Eq. (13. 4p . that 
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may cause quite different statistical weights, as was explicitly shown for the simple example 
in Sec. IIII B[ We can understand this as the appearance of a new, high order, near resonance 
which is no longer disrupted by the bath interaction, resulting in the appreciable occupation 
of the upper (higher average energy) level of that pair, and modifying the statistical weight 
of the states to which it is connected. 

This is emphatically not to suggest the lack of a well-defined statistical behavior of a 
Floquet system with given coupling to a bath. Indeed, the above procedure describes how 
that is to be calculated to any desired accuracy. However, the behavior does depend on 
the specific magnitude and form of that coupling, in sharp contrast to the situation for 
conventional time-independent statistical mechanics. 



V. CONCLUSION 



We have derived the general equations (12.20p for the reduced density matrix, in the 
stationary state, of Floquet systems with weak but finite coupling to a heat bath, previously 
discussed {4, 3, 9, 11 1 only for situations where all quasienergy spacings were large compared 
with the coupling to the bath (and therefore only for systems restricted to a finite number 
of states). New and interesting results give a prescription for practical calculation of the 
steady state statistical properties of many realistic Floquet systems, including ones with 
dense quasienergy spectra. 

If the coupling to the bath F is restricted only to be small compared to the driving 
frequency u and the inverse bath relaxation time I/tc, the diagonal and off-diagonal matrix 
elements are approximately time-independent in the basis of the time-periodic parts of the 
Floquet states. In the special case of a Hilbert space truncated to finite dimensions and 
the coupling weak compared to all quasienergy separations, then the equations simplify and 
are equivalent to those discussed in 11| and Qj. The density matrix is diagonal, but the 
values of the diagonal elements (occupation probabilities of the Floquet states) depend on 
the detailed form (though not on the magnitude, as long as it is sufficiently weak) of the 
coupling of the system to the bath. 

Weakly avoided crossings, where quasienergies become nearly degenerate, are a pervasive 
feature of Floquet systems. Therefore, we have analyzed in detail the general density matrix 
equations in such regions. The rate equations for diagonal density matrix elements in the 
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diabatic basis, Eq. (13.30 . are shown to be identical to those well away from the crossing and 
with arbitrarily weak coupling to the heat bath, except that the avoided crossing introduces 
a new effective rate i?^'^, defined in Eq. (13.41) . between the diabatic states 1 and 2 of the 
crossing. This rate depends on the effective coupling F of the states 1 and 2 to the heat 
bath, the distance d from the avoided crossing, and the relative strength of the coupling to 
the heat bath, F/A. In addition, significant off-diagonal density matrix elements pu and 
P21 may appear (see Eq. (13.61) ). 

The density operator changes drastically with A in the neighborhood of an avoided cross- 
ing, d < 1, if also the bath coupling is weak enough, F ^ A. Quite importantly, significant 
changes occur in the weights not only of the states mixed in the crossing, but in those of 
many other states, as well. At the center of the crossing the two states of the crossing pair 
acquire equal weight, in general a very different relative occupation from that outside the 
crossing region, and this consequently modifies significantly the weights of the network of 
states with which they interact, directly and indirectly, as described by the rate equations. 
This is exhibited explicitly for the simple example worked out in Sec. IIIB, where the occu- 
pations of all states between those labeled rii and 722 are strongly modified in this regime. 
On the other hand, if the coupling to the bath is relatively large, F ^ A, then physical 
intuition suggests that the interaction with the bath effectively broadens the quasienergy 
dispersion curves by more than the gap itself, such that the avoided crossing should have 
no consequences on the statistical properties of the system. We have shown that this is in- 
deed the case, namely that the density matrix (in a A-independent basis) remains virtually 
unaffected. 

Having studied in detail Floquet systems spanned by a finite basis, we have looked at the 
consequences when the size of the basis set increases without limit, as for a particle confined 
by an unbounded potential and driven periodically. We have concluded that, in many cases, 
the important anomalies of such systems, including the appearance of an infinite number of 
weakly avoided crossings and the non- convergence of the Floquet states with increasing basis 
size, play no important role for the physical properties of the system in contact with a heat 
bath, confirming our physical intuition. One can safely employ, in numerical calculations, 
a finite suitably truncated basis set for such a system. However, a special feature of such 
a Floquet system, which has a dense quasienergy spectrum, is that its statistical properties 
depend in general on both the magnitude and the overall strength of the coupling between 
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the system and the bath. 

These calculations have been based on a perfectly periodic driving field. We recognize 
that effects that we have neglected, notably fluctuations in amplitude and frequency of the 
driving fleld, are similarly expected to wash out the effects of sufficiently weakly avoided 
crossings. Those effects will be examined in a future publication. 



Acknowledgement 

We thank S. Kohler and W. Wustmann for helpful discussions. RK thanks the Kavli 
Institute for Theoretical Physics, Santa Barbara for hospitality. DH is grateful for the partial 
support of the Edward A. Dickson Emeriti Professorship and NSF Grant PHY05-51164. 



APPENDIX A: VALIDITY OF LOWEST ORDER BORN APPROXIMATION 



In the still exact Eq. (12. 4p we have made the lowest order Born approximation by using 



the factorization Eq. (12. 5p of the total density matrix directly on the right hand side of that 
equation. In order to check the validity of this approximation, we analyze in this appendix 
the fourth order correction arising from the factorization of the total density matrix only in 
the first iterate of Eq. fl2.4p , namely 
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Each of these two terms has the structure of three nested time integrals, with a sum of 
integrands from the nested commutators (the terms differing only by the order of the factors), 
each of which is a product of 4 system operators and a fourth order correlation function of 
the bath. A representative such bath correlation function is 

{B{t)B{t')B{t")B{t"')) - {B{t)B{t')){B{t")B{t"')) (A2) 
= (^(t)5(t"))(^(t')^(r)) + {mB{t''')){mB{t'')) (A3) 
= G{t - t")G{t' - t'") + G{t - t"')G{t' - t"), (A4) 

where the equality of (A3) to (A2) makes use of the generalized Wick theorem The 
final pair of correlation function products (A4) has the important property that in each of 
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the two terms the time intervals overlap. (The only combination where this is not the case 
is explicitly subtracted in (A2).) Thus all three times appearing in the integrals, t', t", and 
t'", must be close to t in order to give an important contribution, since G{t) is significant 
only for r < Tc. This same condition holds for all contributions to the fourth order term. 
Each such term has an upper limit given by 

T,^j'T!{A'){Br, (A5) 

which compared to the magnitude of the second order rate r2 is small if the weak coupling 
condition, TcT2 -C 1, Eq. fl2.1Up . is fulfilled, namely that during the bath correlation time 
Tc the impact of the second order terms is small; Eq. fl2.10p is itself sufficient to ensure the 
smallness of higher order corrections. 

A ,„o..e .est*«ve cond.ion. P. < s, fo. aU was i,„posed in R.f. Q. na,„e,v *a. the 
system-bath coupling must be smaller than the smallest quasienergy splitting. In contrast, 
the above estimate of the fourth order term shows that larger system-bath couplings with 
r2 > Eij can equally well be treated with the present approach, as long as the weak coupling 
condition, Eq. fl^lB . is fulfilled. 



APPENDIX B: VALIDITY OF TIME AVERAGING OF THE RATES 

The replacement of rates by their averages over a driving period (see Eq. (12.191) and 
the sentence preceding it) relies on the assumption that the density matrix elements do 
not vary substantially over that period, T = 2t^/uo. At first glance this seems to be wrong 
as, e.g., an off-diagonal density matrix Pijit) with an associated quasienergy difference Sij, 
that can be as large as u;, would have in the absence of the coupling to the heat bath a 
substantial variation over the time T. But it will turn out that such off-diagonal matrix 
elements (corresponding to large quasienergy differences) are in fact small at large times. 

We use the observation that the set of differential equations (I2.18P is linearly coupled 
and has time-periodic coefficients Rij-ki{t) with period T. Thus one can apply the Floquet 
theorem, which requires the general solution to be of the form 

n 

with time-periodic functions rij.n{t + T) = rij.n{t), and e'*'""^ an eigenvalue of the time- 
evolution operator for all density matrix elements over one period. (The A„ are assumed 
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here to be non-degenerate.) This solution is in general not periodic in time. 

In the following we are interested in the solution for large times. All solutions with 
Re{Xn) < decay to zero in this limit and are of no further interest. From the general 
properties of a density matrix, < p^^ < 1 and |pijp < PaPjj < 1, it follows that solutions 
with Re{\„) > are impossible. So we are left with discussing solutions with _Re(A„) = 0. 
For the trace of the density matrix to be 1, i.e. time-independent and nonzero, there must 
be a solution Aq = 0. Thus we have a stationary solution for the density matrix elements 

Pij{t)=ri,.o{t), (B2) 

which is time-periodic, with at least some diagonal components which have non-zero time 
averages. (We cannot exclude that in addition to this, there is a traceless quasiperiodic 
solution from pairs of purely imaginary A„. As its time average would be zero, we ignore it 
in the following.) 

We now expand this stationary time-periodic solution in Fourier coefficients 

K 

SO that the above differential equations ( I2.18P reduce, for each of the Fourier components, 
to 

-i{e,j+Kuj)p,,{K) = [pijiK - M)R,k-ik{M) + Pii{K - M)R*,.i,{-M) (B4) 

^ k,l,M 

- pu{K - M){Ri,,,,{M) + Rl,,i^{-M)) 

The right hand side is second order in the small system-bath coupling constant 7, because 
each rate factor fl2.15p is. Although we are explicitly concerned with couplings which are 
not arbitrarily weak with respect to the smallest quasienergy differences in the problem, we 
do want to consider, as in the conventional static case, coupling which is weak on higher 
energy scales, such as the driving frequency uo. Then for non-zero values of K the necessary 
smallness of the left hand side of flB4l) must be achieved by the smallness of pij{K), and it is 
consistent to neglect all such components of the density matrix in the solution; we will set 

p,,(ir^O) = 0, (B5) 

and we will simply write pij for Pij{K = 0), which is by definition time independent. (We 
point out that this implies a time-periodic density operator for large times, as one would 
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expect physically.) This time independence of pij at large times then implies in (lB4p that 
only M = terms appear — i.e., that the complex rates R can be replaced by their time 
averages R{M = 0). (Note that for a quasienergy difference eij ~ ±uj between a Floquet 
state close to the lower boundary and a state close to the upper boundary of the chosen 
strip < e < uj, one would need to choose K = =pl in Eq. (lB4p or choose a conveniently 
displaced strip.) 

APPENDIX C: DERIVATION OF THE EFFECTIVE RATE i?^^ 

In this appendix we derive the rate equations for the reduced density matrix elements 
in the diabatic basis of an isolated weakly avoided crossing in the stationary limit. This 
leads to an effective rate R^'^, which expresses the entire A-dependence related to an avoided 
crossing. 

Starting from the operator equation (12. 8p one follows the steps leading to Eq. (12.201) . but 
now in the diabatic basis of the avoided crossing, which is related to the adiabatic basis by 

\ua{t)) = a\u,{t)) + (3\u2{t)) (CI) 
M)) = I3\u,{t)) - a\u2it)) (C2) 

with a = \J (1 + (i/Vl + c?^)/2, 13 = — + d?)Ji and d, the dimensionless distance 

from the center of the avoided crossing, is defined in Eq. (13. 2p . This basis change alters the 
left hand side of Eq. (I2.20p . as the coherent dynamics is less conveniently described in the 
diabatic basis. The right hand side is left invariant, but now with rates in the diabatic basis, 
and can be simplified with the following assumptions: 

1. We assume that the avoided crossing is sufficiently weak and isolated within the 
quasienergy spectrum, such that we can neglect any dependence of diabatic states 
1 and 2 on the driving amplitude A, valid over a sufficiently limited range of A near 
Aq. We assume the same for all other Floquet states, such that the diabatic basis is 
independent of A within that range. 

2. All other quasienergy separations near A = Aq are assumed to be much larger than 
the splitting A as well as larger than the effective coupling to the bath, such that all 
off-diagonal density matrix elements can be neglected, except for pi2 and p2i- 
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3. We define rates Rij-ki in the diabatic basis in analogy to Eqs. (12.190 and (12.151) . where 
the A dependence appearing in the quasienergies ei and 62 is neglected in the neigh- 
borhood of the avoided crossing. That is, the Fourier transformed bath correlation 
function g{s — muj) is taken to be constant over this very small range of e. This is a 
good approximation as long as Ar^ ^ 1 holds, which is the case for the weakly avoided 
crossings we are interested in. Together with assumption 1 this makes the rates Rij-^ki 
independent of A. 

4. We set Ai2{m) = A2i{m) = (see Eq. (12.131) ). as we assume that the spatial structures 
of states 1 and 2 are very different. We are explicitly looking at a very weakly avoided 
crossing, and this typically arises from such dissimilar pairs of states, such as those 
arising from undriven states one of which is of low and the other of very high energy. 
Similarly, any other state j cannot be spatially similar to 1 and 2 simultaneously, 
leading to Aij{m)A2j{m) = 0. 

5. For simplicity we restrict the discussion to periodic drivings that are symmetric, i.e. 
Hs{t) = Hs{~t). This makes the rate -R22;ii; which appears in the definition of F in 
Eq. ([33]), real. 

For the density matrix in the subspace of states 1 and 2 we find 

A Im(pi2) = -pu^Rik + ^pkkRki (C3) 

k k 

-A Im(pi2) = -p22^R2k + ^PkkRk2 (C4) 



iA 



(pii - P22) - 2(ipi2 = pi2r (C5) 



with the symbol F introduced in Eq. (13. 5p . 

The third equation expresses the off-diagonal element pi2 in terms of the diagonal ele- 
ments, as explicitly written in Eq. (13. 6p . Consequently, the left hand side of the first two 
equations can be written in terms of the diagonal elements, 

Alm(pi2)=i?"^(pii-P22) (C6) 



where we use the notation R^^ introduced in Eq. (13.41) . We now redefine the effective rates 
Rij in the diabatic basis, equal to the original rates Rij except for R12 = R^'^ and similarly 
R21 = R^^ (note R12 = by assumption 4 above). This allows us to rewrite the first two 
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equations in the standard rate equation form, as if there were just diagonal density matrix 
elements: Eq. ( \3.3\i . 



APPENDIX D: DIAGONALIZING THE DENSITY MATRIX 



It is instructive to ask in which basis the reduced density matrix is diagonal near a weakly 
avoided crossing, and whether this is close to the diabatic 1,2 basis above, or close to the 
a, b basis of the Floquet states. The relevant sector of the density matrix in the 1, 2 basis is 
given by 

Pii A{Pn - P22) 
^^*(Pii-P22) P22 

with 

A 



P 



1 



2d-iV//\' 



A 
\A\ 



It is diagonal in the basis of the two states 
l^di) = 

\Ud2) = 

with 

leading to diagonal elements 



1 + - \ui) + Jl-- \U2) 

z \ z 



/V2 

/V2, 



1 + 4|A|2, 



(Dl) 



(D2) 



(D3) 



(D4) 

(D5) 
(D6) 



Pdl,dl = [pil(l + 2;) + P22(l - ^)]/2 
Pd2,d2 = [Pll(l - 2;) + P22(l + ^)]/2. 

One can check that this gives the expected result in two limiting situations: 

• Far from the avoided crossing, d ^ 1, 01 for F ^ A , the effective rate R^'^ between 
states 1 and 2 becomes small, we have A —>■ and z — > 1, and we regain the 1,2 basis 
as that in which the reduced density matrix is diagonal. 

• For weak coupling, F -C A and near the avoided crossing, c? <^ 1, we have z ^ 1. It is 
immediately clear that at the center, where d ^ 0, the appropriate basis states for a 
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diagonal density matrix become nearly the equal admixtures of states 1 and 2, which 
are the Floquet states a and h. In fact, this is the case for all d in this limit of weak 
coupling to the bath; for infinitesimal F these states which diagonalize p are exactly 
the Floquet states. 

We can readily understand these two limits in terms of the balance between energy 
pumping of the system by the periodic field and relaxation by the bath. For F <g A the 
near resonant periodic driving dominates, and the Floquet states in the absence of the bath 
are approximate solutions, only weakly relaxed. In the opposite limit, F ^ A, those Floquet 
levels are effectively broadened more by the bath interaction than their separation, and the 
mixing of the diabatic states is a perturbation on their relaxation by the bath. 
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